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1. INTRODUCTION 

Reflector antennas on earth orbiting spacecrafts generally 

cannot be described analytically. The reflector surface Is 

subjected to large temperature fluctuation and gradients, and 

Is thus warped from Its true geometrical shape. Aside from 

distortion by thermal stresses, reflector surfaces are often 

purposely shaped to minimize phase aberatlons and scanning 

losses. To analyze distorted reflector antennas defined by 

discrete surface points, a numerical technique must be applied 

to compute an Interpolator surface passing through a grid of 

discrete points. Although numerical techniques for analyzing 

reflector surface distortions and doubly curved reflector sur- 
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faces have been reported In open literatures, all these 
techniques are rather complicated and Involve lengthy computa- 
tions. In this paper, a simple numerical technique based on 
Taylor's expansion and finite differences approximation Is 
presented. By applying the numerical technique to approximate 
the surface normals of a distorted reflector surface, the aper- 
ture field over a near-field plane Is obtained using geometric 


optics technique. The aperture Integration method Is then 
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applied to find the radiation patterns. 

2. FORMULATION 

The proposed numerical technique Is a quadratic approxlma- 
tlon of a distorted surface defined by a grid of nine points, 
Let (x,y,z) be the coordinate of any point, P, on the surface 
which Is described by 

z * z(x,y) 

The position vector of P Is given by 


r = XX + yy + zz 


To obtain an approximation for the surface z(x,y), the position 
vector Is replaced by Its Taylor's expansion around r (x Q , 
y o * z q) * c ^ oses ^ P°^ n t to the point of Incidence. That Is 


t(i.,y,z) = r 0 (x 0 ,y 0 ,z 0 ) + <* - *„> fj ♦ <y - y„) ay 
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where the partial derivatives of r are evaluated at r Q (x 0 , 

y 0 ,z Q ). All terms higher thai> the quadratics of x and y 
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have been neglected as Indicated by 0(x ,y ). The surface 
may then be approximated by 


z(x,y) » z • t « z Q + (x - x 0 ) — + (y 
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All the partial derivatives of z are evaluated at "^ 0 ( x Q , 
y 0 ,z 0 ). The normals to the Interpolatory surface can be 
found by taking the gradient of z(x,y) as follows: 


"x ■ 8 * < x - x o> J ♦ < y - V Tk * 0(x2 ’ y2) 

"y ’ If * <x - x o> * < y ' V *• 0<x2 ' y2) 


n z ■ - 1 


where (n ,n ,n ) are the x, y, and z components of n with 

a y l 

the normal pointing toward the source. Finite differences can 
be applied to approximate the partial derivatives. Since the 
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partial derivatives are of second order, a grid of nine points 
Is required. 

The partial derivatives may be expressed as 



< z 4 - z 0> < z 0 - z 5> 

<*4 - X 0> ~ ftp - X S> 

( X 4 ~ X S ^ 

2 


3Z 

ax 



z 

X 


5 

5 



(Z, 


Z 0 ) 


<Z r 


Z 7> 


(y 2 - y 0 ) (Vp - y 7 ) 

( yo ~ 


32 Z 2 ~ Z 7 

ay ~ y 2 - y 7 

^ Z 1 ~ z 3 ^ ^ Z 6 ~ Z 8 ^ 
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where (x 2 ,y 2 ,z 2 ) ... ( x 8 *y 8 z 8 ) are the coordi- 

nates of the points around the closest point ( x Q , y Q , z Q ) - 

3. DISCUSSION AND RESULT 

In applying the numerical technique to compute the surface 
normals, the nine closest points to the point of Incidence are 
located by searching an array of discrete surface points defin- 
ing the distorted antenna surface. Ray-tracing Is then used to 
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compute the reflected field from the approximated normals. The 
secondary far field patterns are obtained from the Fourier 
transformation of the aperture field distribution. 

To Illustrate the numerical technique, the numerical algo- 
rithm has been programmed and applied to the antenna configura- 
tion shown In figure 1. The surface grid points for an 
artificially "enlarged" parabolic surface are generated. Thus, 
the boi'i’sdary points along the rim of the actual reflector sur- 
face can be Included In the computation. The computed secondary 
far field patterns for the symmetric and the offset cases are 
shown In figures 2 and 3 respectively. In the figures, the 
E-plane radiation patterns computed from an analytic expression 
of a parabolic surface are drawn with solid lines, while that 
computed numerically are drawn with dots and broken lines. The 
E-plane radiation patterns computed numerically from a grid of 
50x50 surface points of ~0.7 x spacing was found to compare 
exactly to that computed from an analytic expression. For the 
offset case, the aperture surface Is described by fewer points; 
therefore, a higher density of target points Is required to 
produce the same accuracy. With a scan step of \/3 for the 
symmetric case and >/6 for the offset case, the computer time 
required to compute the secondary radiation patterns numerically 
Is 25 sec for the symmetric case and 17 sec for the offset case, 
respectively. For comparison purpose, the radiation pattern 
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for the offset case has been computed using spline function 

5 

approximation. The computed E-plane pattern Is shown In 
figure 3. 
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